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Abstract 

We study the moments of hadronic event shapes in e + e _ annihilation within the context of next-to- 
next-to-leading order (NNLO) perturbative QCD predictions combined with non-perturbative power 
corrections in the dispersive model. This model is extended to match upon the NNLO perturbative 
prediction. The resulting theoretical expression has been compared to experimental data from JADE 
and OPAL, and a new value for a 3 (Mz) has been determined, as well as of the average coupling ao 
in the non-perturbative region below [ii — 2 GeV within the dispersive model: 

a s (M z ) = 0.1153 ± 0.0017(exp) ± 0.0023(th), 
a = 0.5132 ± 0.0115(exp) ± 0.0381(th) , 

The precision of the a s (Mz) value has been improved in comparison to the previously available next- 
to-leading order analysis. We observe that the resulting power corrections are considerably larger 
than those estimated from hadronization models in multi-purpose event generator programs. 



1 Introduction 



Event shape variables measure geometrical properties of hadronic final states at high energy particle 
collisions. They have been studied extensively at e + e~ collider experiments, which provided a wealth 
of data at a variety of centre-of-mass energies. Exploiting this large energy range, one can attempt 
to disentangle perturbative and non-perturbative contributions (which scale differently with increasing 
energy) to event shape observables. 

Apart from distributions of these observables, one can also study mean values and higher moments. 
The nth moment of an event shape observable y is defined by 



where y max is the kinematically allowed upper limit of the observable. Moments were measured for a 
variety of different event shape variables in the past. The most common observables y of three-jet type 
are: thrust T [1] (where moments of y = (1 — T) are taken), the heavy jet mass p — M\js [2], the 
C-parameter [3] , the wide and total jet broadenings B w and B T [4] , and the three-to- two-jet transition 
parameter in the Durham algorithm Y 3 [5]. Definitions for all observables are given in, for example, 
Ref. [6]. Moments with n > 1 have been measured by several experiments, most extensively by JADE [7,8] 
and OPAL [9], but also by DELPHI [10] and L3 [11]. A combined analysis of JADE and OPAL results 
has been performed in Ref. [12]. 

As the calculation of moments involves an integration over the full phase space, they offer a way of 
comparing to data which is complementary to the use of distributions, where in general cuts on certain 
kinematic regions sxe c ipplied. Furthermore, the two extreme kinematic limits two-jct-like events and 
multi-jet-like events - enter with different weights in each moment: the higher the order n of the moment, 
the more it becomes sensitive to the multi-jet region. Therefore it is particularly interesting to study 
the NNLO corrections to higher moments of event shapes, as these corrections should offer a better 
description of the multi-jet region due to the inclusion of additional radiation at parton level. 

Moments are particularly attractive in view of studying non-perturbative hadronization corrections 
to event shapes. In event shape distributions, one typically corrects for hadronization effects by using 
generic Monte Carlo event simulation programs. A recent study, carried out in the context of a pre- 
cision determination of the strong coupling constant from event shape distributions [13], revealed large 
discrepancies between the standard event simulation programs used at LEP [14, 15] on one hand and 
more modern generators [16], which incorporate recent theoretical advances, on the other hand. In the 
event shape distributions, it is very difficult to disentangle hadronization corrections empirically, since 
they typically result in a distortion of the distribution, which can not be unfolded in a straightforward 
manner. 

In event shape moments, one expects the hadronization corrections to be additive, such that they can 
be divided into a perturbative and a non-perturbative contribution, 



where the non-perturbative contribution accounts for hadronization effects. Based upon the calculation of 
next-to-next-to-leading order (NNLO) QCD corrections to the event shape distributions, which became 
available recently [6, 17-21], the perturbative contribution to event shape moments is now known to 
NNLO [22,23]. The non-perturbative part is suppressed by powers of X p /Q p (p > 1), where Q = \fs is 
the centre of mass energy and Ai is of the order of Aqcd- The functional form of X p has been discussed 
quite extensively in the literature, but as this parameter is closely linked to non-perturbative effects, it 
cannot be fully derived from first principles. 

In this work, we use the dispersive model derived in Ref. [24-27] to compute hadronization correc- 
tions to event shape moments. This model provides analytical predictions for the power corrections, 




(1) 



(y n ) = (y n ) P t + (y n ) 



(2) 
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and introduces only a single new parameter eto, which can be interpreted as the average strong coupling 
in the non-perturbative region. This model has been used extensively in combination with NLO QCD 
perturbative calculations to study event shape moments [9,28-30]. To combine the dispersive model with 
the perturbative prediction at NNLO QCD, we extended its analytical expressions to compensate for all 
scale-dependent terms at this order. By comparing the newly derived expressions with experimental data 
on event shape moments, we perform a combined determination of the perturbative strong coupling con- 
stant a s and the non-perturbative parameter ao- Compared to previous results at NLO, we observe that 
inclusion of NNLO effects results in a considerably improved consistency in the parameters determined 
from different shape variables, and in a substantial reduction of the error on a s . 

In Section 2, wc outline the structure of perturbative and non-perturbative contributions to event 
shape moments. The predictions of the dispersive model to power corrections are extended to NNLO 
in Section 3, and used to extract a s and ao from experimental data in Section 4. In Section 5 the 
results obtained within the dispersive model are compared to those from multi-purpose event generator 
programs. 

2 Power corrections to event shape moments 

Non-perturbative power corrections can be related to infrared rcnormalons in the perturbative QCD 
expansion for the event shape variable [24,25,31-36]. The analysis of infrared renormalon ambiguities 
suggests power corrections of the form X p /Q p , but cannot make unique predictions for X p : it is only the 
sum of perturbative and non-perturbative contributions in (2) that becomes well-defined [37]. Different 
ways to regularise the IR renormalon singularities have been worked out in the literature [38-43] . 

One approach is to introduce an IR cutoff p,i and to replace the strong coupling constant below the 
scale by an effective coupling such that the integral of the coupling below p,j has a finite value [24-27] 

— / dQa cS {Q 2 ) =a {m) . (3) 
Mi Jo 

This dispersive model for the strong coupling leads to a shift in the distributions 

where the numerical factor a y depends on the event shape and is listed in Table 1, while P is believed 
to be universal (universality breaking terms arise from hadron mass effects [44] in the moments of p, 
an estimate on these effects can be obtained from general-purpose event generator programs, e.g. from 
PYTHIA [14]) and scales with the CMS energy like m/Q. 

By inserting (4) into the definition of the moments, one obtains: 

/■2/max 1 da 

(V n ) = / dyy n — (y) (5) 
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Table 1: The a y coefficients of the non-perturbative event shape moment prediction 
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discarding the integration over the kinematically forbidden values of y. This leads to the the non- 
perturbative predictions for the moments of y: 

(2/ 2 > = (y 2 )pt + 2(y 1 )p t KP) + (a !y P) 2 , 

(y 3 ) = (y 3 ) P t + 3(y 2 ) pt (a y P) + 3(yX t (a v P) 2 + (a v P) 3 , 

(y 4 ) = (y 4 ) P t + 4(j/ 3 ) pt KP) + 6(y 2 ) pt (a y P) 2 + 4(j / 1 ) pt KP) 3 + {a y P)\ 

(y 5 ) = (y\t + H(y\t(a y P) + lO(y\ t (a y P) 2 + 10(y 2 ) pt (a y P) 3 + ^yX^Pf + (a y Pf (8) 

It should be noted that the multiplicative power correction in (7) is considered to be accurate to l/Q- 
For n > 2, the evaluation (8) yields also higher powers of P, which are formally of higher order in inverse 
powers of Q. Contributions with the same scaling behaviour could equally come from subleading power 
corrections in P. Compared to the terms above, these subleading power corrections would be weighted 
with higher perturbative moments, and are thus suppressed numerically. 

The perturbative contribution to (y n ) is given up to NNLO in terms of the dimensionless coefficients 
*A.y,m $y,n and C y _ n as. 

<y">pt(^ 2 ) = (^)A.«+(^) 2 (^ + A,^oiog^ 

^) 3 (c y , n + 2B y , n /3 log ^ + A y , n (p 2 log 2 — log ^ 

Z7T J \ S \ S S 

+0(ai) . (9) 

In here, s denotes the centre-of-mass energy squared and n is the QCD renormalisation scale. The 
NLO expression is obtained by suppressing all terms at order o? s . The first two coefficients of the QCD 
/3-function are 

\\C A - AT R N F 
6 

17C 2 - IQC A T R N F - 6C F T R N F 
Pi = g , (W) 

with C A =N,C f = (N 2 - 1)/(2N), T R = 1/2 for N = 3 colours and N F quark flavours. 

The perturbative coefficients in (9) are independent on the centre-of-mass energy. They are obtained 
by integrating parton-level distributions, which were calculated recently to NNLO accuracy [6,17,19]. 
These parton-level calculations are based on a numerical integration of the relevant three-parton, four- 
parton and five-parton matrix elements, which are combined into a parton-level event generator [18,20,21] 
after subtraction of infrared singular configurations using the antenna subtraction method [45] . These 
NNLO event shape distributions were used subsequently for improved extractions of the strong coupling 
constant [13,46-49], matched on all-order resummation of logarithmically enhanced corrections [48,50], 
and used for power correction studies on the thrust distribution [47]. 

The coefficients entering the event shape moments are computed at a renormalisation scale fixed to the 
centre-of-mass energy, and are therefore just dimensionless numbers for each observable and each value 
of n. For the first five moments of the six event shape variables considered here, they were computed up 
to NNLO in [22,23]. 



3 Dispersive model extended to NNLO 

Up to now, the dispersive model for power corrections to event shapes was used in connection with NLO 
calculations of the perturbative part. In this context, one obtains the following, l/Q-dcpcndcnt power 
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correction [26]: 



P = ■ M ■ i a - 



a s { m ) + ^ al( m ) fln^ + 1 + A ) + 0{p? B ) 
7T \ Mr 2 Po 



with the Milan factor = 1.49 ± 20%, which is known at two loops. Its uncertainty [51] accounts 
for currently unknown corrections beyond this loop order. The term in square brackets amounts to the 
renormalon subtraction in the power corrections, expanded to NLO. 

The prediction of the dispersive model can be extended to match onto the NNLO perturbative pre- 
diction, and first steps in this direction were taken already in [47] for power corrections to the thrust 
distribution. 

The perturbative ingredients to the dispersive model are the running of the coupling constant and 
the relation between the MS-coupling and the effective coupling, whose definition [52] absorbs universal 
correction terms from the cusp anomalous dimension. 

In the present context, we use the evolution of the coupling constant to two loops 
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Moreover, the relation between MS-coupling and effective coupling reads 
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(13) 
(14) 

(15) 



The coefficient L is obtained from the three- loop cusp anomalous dimension [53,54], which can be ex- 
tracted from the three-loop corrections to the partonic splitting functions [55] or to the quark and gluon 
form factors [56]. 

The derivation of a generic power correction starts from considering a dimensionless quantity 



F 



(16) 



with 



f(fj,) oc a F a a (n) 



(17) 



assuming F to be dimensionless. The value of p determines the scaling behaviour of the power correction, 
with p = for the leading power correction to event shape variables. 

The dispersive model assumes that in the non-perturbative range of (16) the perturbative strong 
coupling a s (/i) is replaced by an effective coupling that remains finite for all /x values. One defines then 
the value of the integral over this region by 



p 



TtMmt) 



(18) 
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introducing an infrared matching scale ^j,Aqcd <C /«/ <C Q and a p as a non-perturbative parameter. 
One has then to subtract the perturbative part of (16) in the range from to /i/ from the whole integral, 
that is, the value of (18) with a Sj iR replaced by a s . 

This perturbative contribution to (16) thus acquires a dependence on the renormalisation scale /ir used 
in the strong coupling constant. By requiring F to be scale- independent, one can then infer logarithmic 
terms in the non-perturbative contribution to (16). Applied to the event-shape power correction P (with 
p = 0), this results in 



P = 
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(19) 



Together with (9) this gives the full expression for the event shape observable moments, including per- 
turbative and non-perturbative contributions. 

For Bt and B\y there is a further correction to (19). It arises from the kinematical mismatch between 
parton direction and thrust direction used to define the hemispheres used in the broadening variables. 
Retaining (8), this modification can be accounted for by a modification to the power correction. In [27], 
this modification was computed to NLO for the first moment as 



P(B W ) = P 



8C F a s (l + ^ 



3 
4 



A. 

6C F 



Vo 



(20) 



P (Bt) = P 
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with a s (Q) = a s (e~*Q) and 770 = —0.6137. Corrections to higher moments have not been derived up 
to now, and we assume that they can be approximated by using the above modifications to the power 
correction in all moments. The full NNLO expression for these has not been calculated either. The 
potentially dominant NNLO terms can however be approximated by including the effective coupling to 
this order, resulting in 



(Bw) 
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(23) 



However, further NNLO corrections to this expression will reside in the coefficient 770. Therefore, we will 
treat Bw and Bt separately from the other variables in the numerical studies in the following section. 
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nominal value 


up variation 


down variation 


MGcV] 


2 


3 


1 




1 


2 


0.5 


M 


1.49 


1.788 (+20%) 


1.192 (-20%) 



Table 2: Table of the pj, x^ and M variations 



4 Analysis of JADE and OPAL data 

The theoretical expressions for event shapes derived in the previous section contain two parameters: the 
strong coupling constant a s (Mz) and the non-perturbative coupling parameter ao- Using experimental 
data on event shape moments, it is possible to fit these parameters. The data from the JADE and OPAL 
experiments [8] consists of 18 points at centre-of-mass energies between 14.0 and 206.6 GeV for the first 
five moments of T, C, Y-j, My, B\y and Bt, and have been taken from [29]. For each moment the NLO 
as well as the NNLO prediction was fitted with a s (Mz) and ao as fit parameters, except for the moments 
of Y3, which have no leading ^ power correction and thus are independent of ao. For the heavy jet mass, 
we use only the even moments (M^) and (Mjj-), since the theoretical prediction is in terms of p = M\js. 

4.1 Fits 

The fits were done using the program ROOT [57] and its x 2 fit method. The errors used for the fit 
were the total errors, composed of the experimental statistic and systematic errors, added in quadrature. 
Based on these, ROOT returned errors on the fit which are displayed in Tables 4-15 in the appendix 
together with the fit results. For T and C the NNLO values of a s (Mz) and a seem to be more stable 
throughout the moments, as at NNLO they increase less towards higher moments than at NLO. For Y3 
and p, where the values decrease at higher moments, this is not the case. These moments show a s {Mz) 
results which are significantly lower at NNLO than at NLO. The a s (Mz) values of B\y are much lower 
than the ones of the other observables, and do not change much from NLO to NNLO. For Bt the a s (Mz) 
values at NNLO are lower than at NLO. Both are exceptionally stable throughout the different moments. 
The ao values of all moments are higher at NNLO than at NLO. 

4.2 Theoretical systematic errors 

There are different parameters in the theoretical prediction which may influence the results displayed 
above, namely the matching scale pi, the renormalisation scale pr and the Milan factor M.. In order to 
estimate the resulting theoretical uncertainty on a s (Mz) and ao, the fits were repeated, pi, pn and M. 
being separately varied by a certain amount. 

For this purpose the scaling factor x^ = ^ was introduced. The uncertainty on the corresponding 
parameter was then taken to be the difference between the nominal and the new value returned by 
ROOT. In order to get a total systematic error, the greater values of the up and down uncertainties were 
determined and quadratically added. As ao depends directly on pj no error was determined for this 
variation. For Y 3 there is only an error on a s (Mz) coming from the variation, since the theoretical 
description of this observable does not contain a contribution from the leading power correction, and is 
thus independent on pj and A4. At NLO the fit to the moment (C 3 ) suffers from a numerical instability 
by scaling up M. by 20%. The numbers reported in Table 6 refer to an up variation of 19%. 

The NLO error on a s (Mz) agrees well with the values of [8]. At NNLO, it is reduced by more than 
half throughout all event shape observables except Byy, confirming a good description by the NNLO 
prediction. Unfortunately, this is not the case for the error on ao- It does not change much from NLO 
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Figure 1: Plot of the individual measurements for a s (Mz). The shaded region corresponds to the error 
band defined by the weighted mean for the event shape and the total error on it. 



to NNLO, even increasing a little in the first moments due to the higher x^ uncertainty at NNLO and 
decreasing slightly at the higher moments, with exception, again, of Bw- Analysing the different sources 
of the systematical errors, we observe that the error on a s (Mz) is clearly dominated by the x^ variation, 
while the largest contribution to the error on ckq comes from the uncertainty on the Milan factor M . Since 
this uncertainty has not been improved in the current study, it is understandable that the systematic 
error on ao remains unchanged. This finding clearly motivates the need for a three-loop calculation of 
the Milan factor. However, it is very important to note that the uncertainty on the Milan factor has little 
impact on the extraction of a s (Mz), thereby demonstrating the systematic decoupling of perturbative 
and non-perturbative effects in the dispersive model. 

For the higher moments (n > 2) of the jet broadenings Bw and Bt, the kinematical modifications 
to the power correction are not known at present. We have approximated them in the above fits by 
the corrections to the first moments, given to NLO and NNLO in the previous section. If we do not 
apply these correction to the higher moments, the mutual consistency of the parameter extractions from 
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Figure 2: Error band plot of the individual measurements for ao 



different moments of Bt deteriorates considerably, while only minor improvements in consistency are 
observed on Bw 

Including empirical hadron mass corrections [44,58] from PYTHIA affects in particular the parameter 
extraction from p, resulting in values of a s (Mz) and ao from p much lower than from the other variables. 
Since these corrections may interplay with other non-perturbative parameters in PYTHIA, we do not 
include them in our default fits or error estimates. 

By taking the weighted means over the corresponding values from all moments of all observables 
one gets combined values for a s (Mz) and ao- The weights are given by the inverse of the total error 
squared and are normalized such that the sum over all weights is equal to one. For the errors one has 
to take care of the correlation between the errors of the single measurements. The correlation matrix for 
a s (Mz) and ao is in first approximation equal to the correlation matrix for the event shape moments, 
since the variable transformation is linear in first approximation. The correlation matrix for the event 
shape moments is given in [29] . We first combine the measurements from different moments of the same 
observable. Figures 1 and 2 compare the combined NNLO results on the a s (Mz) and ao measurements. 
Owing to the large correlation between individual moments of the same observable, the combined errors 
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Theoretical Error 
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Experimental Error 
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r 


0.5444 


0.0184 


0.0388 


0.0430 


C 


0.4841 


0.0066 


0.0347 


0.0353 


P 


0.6380 


0.0270 


0.0824 


0.0867 


B T 


0.4924 


0.0102 


0.0449 


0.0460 


Bw 


0.3362 


0.0125 


0.0338 


0.0360 


Total 


0.4604 


0.0108 


0.0359 


0.0375 


Total w/o Bt,Bw 


0.5132 


0.0115 


0.0381 


0.0398 



Table 3: Table of the a s (Mz) and «o results for the individual moments and the global weighted average. 

are only marginally smaller than the errors obtained from single measurements. The combined results 
and their errors are summarised in Table 3. From this Table, we clearly observe that the the theoretical 
error on the extraction of as(Mz) from p, Y$ and B\y is considerably smaller than from r, C and Bt- 
It was observed previously in [22] that the moments of the former three shape variables receive moderate 
NNLO corrections for all n, while the NNLO corrections for the latter three are large already for n = 1 
and increase with n. Consequently, the theoretical description of the moments of p, Y-j and Bw displays 
a higher perturbative stability, which is reflected in the theoretical uncertainty on a${Mz) derived from 
them. 

In a second step, we combine the a s (Mz) and ag measurements obtained from different event shape 
variables. Taking the weighted mean over all values, but excluding the values for the moments of Bw 
and Bt where the theoretical description is incomplete, we obtain at NNLO: 

a s (M z ) = 0.1153 ± 0.0017(exp) ± 0.0023(th), 

a = 0.5132 ± 0.0115(exp) ± 0.0381(th) , (24) 

where the errors have been derived taking into account the correlation between the moments of different 
event shapes. Including the values for Bw and Bt modifies this result to: 

af(M z ) = 0.1131 ± 0.0017(cxp) ± 0.0022(th), 
= 0.4604 ± 0.0108(exp) ± 0.0359(th) . 

These latter values are however quoted only to illustrate the impact of including the broadenings. The 
default fit result is (24), where only observables with a consistent theoretical description are included. 

To illustrate the improvement due to the inclusion of the NNLO corrections, we also quote the 
corresponding NLO results. Based on r, C, p and Y3, we obtain: 

a™ LO (M z ) = 0.1200 ± 0.0021(exp) ± 0.0062(th), 
c^ LO = 0.4957 ± 0.0118(exp) ± 0.0393(th) , 
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while inclusion of B\y and Bt modifies this to 



a™' B (M z ) = O.f 147 ± 0.0020(exp) ± 0.0046(th), 
a NLO,s = 40ig ± .0130(exp) ± 0.0296(th) , 

We compare the NLO and NNLO combinations in Figure 3. It can be seen very clearly that the mea- 
surements obtained from the different variables are consistent with each other within errors. The average 
of a s (Mz) is dominated by the measurements based on p and I3, which have the smallest theoretical 
uncertainties. From NLO to NNLO, the error on a s (Mz) is reduced by a factor two, and the result shifts 
towards the lower end of the NLO error band, as was already the case in the individual measurements. 
No improvement and no shift in the central value between NLO and NNLO is seen on a . 

5 Comparison with PYTHIA hadronization corrections 

The primary motivation for studying power corrections to moments of event shapes in the dispersive 
model comes from the observation that the commonly used method to derive hadronization corrections 
from multi-purpose event generator programs may be unreliable [13]. 

To quantify the difference of both approaches to hadronization corrections, we compare them on the 
example of the moments of 1 — T. For this comparison, we extracted the PYTHIA [14] hadronization 
corrections to these moments from the ratio of PYTHIA hadron level and parton level results. Using 
these corrections in combination with the NNLO perturbative expressions for the event shape moments, 
we repeated the fit of as(Mz) on the different moments of 1 — T. The results are displayed and compared 
with the fits in the dispersive model in Figure 4. We observe that both approaches yield a reasonable 
description of the experimental data, but that the resulting values of a s {Mz) are considerably larger when 
applying hadronization corrections extracted from PYTHIA. Given that the perturbative contribution 
increases monotonously with a s (Mz), this indicates that the hadronization corrections in PYTHIA are 
considerably smaller (and perhaps underestimated) than those obtained in the dispersive model. 

This observation is quantified on the moments of (1 — T) at y/s = Mz, displayed in Figure 5. 
Depending on the moment number, we observe that the PYTHIA hadronization corrections are between 
two and four times smaller than those obtained from the dispersive model. It can also be seen that 
the PYTHIA-based predictions are systematically below the experimental data, which perhaps indicates 
that, despite the decent agreement on the full range of energies, Figure 4, PYTHIA fails in the precise 
description of the energy dependence of the hadronization corrections. 

Our comparison suggests strongly that hadronization corrections extracted from PYTHIA (or from 
other comparable multi-purpose event generator programs [15]) are lower than power corrections obtained 
from analytical hadronization models. This can be partially understood from the fact that perturbative 
predictions of PYTHIA are larger than fixed order calculations since they are based on a parton-showcr 
approximation to the all order resummation result, which is known to shift the distribution away from the 
two jet singularity. As a consequence, using PYTHIA hadronization corrections to analyse data in view 
of precision extractions of a s (Mz) may result in anomalously large values, since the missing numerical 
magnitude of the power corrections must be compensated by a larger perturbative contribution. 

6 Conclusions 

In this paper, we studied the perturbative and non-perturbative contributions to the moments of event 
shapes in e + e~ annihilation. In view of the recently calculated NNLO perturbative contributions [22,23] 
to the event shape moments, we extended the dispersive model for non-perturbative power corrections [24- 
26] to include all logarithmic corrections to this order. The normalisation of the power correction (the 
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Figure 3: Error band plot of the final results. The points for a 8 (Mz) are C, T, Y3, Mh and for ao C, T, 
M H . 

Milan factor [26]) is however still restricted to NLO accuracy, and specific corrections [27] to the jet 
broadenings Bw and Bt are also only included to NLO. 

We used this newly obtained theoretical description of the event shape moments to reanalyse data 
from JADE and OPAL in view of a determination of the strong coupling constant ct s {Mz) and of the non- 
perturbative parameter ao . We observed that inclusion of the NNLO corrections results in a considerably 
better consistency among the values extracted from different moments of the same variable, and an 
improved consistency among the different variables. Averaging over the different moments and different 
shapes (excluding By/ and B?, where the theoretical description is incomplete, and taking proper account 



11 






Figure 4: Comparison of fits with hadronization corrections from PYTHIA and power corrections from 
the dispersive model. 



of the uncertainty due to missing terms in the Milan factor), we obtain the following combined values: 

a s (M z ) = 0.1153 ± 0.0017(fit) ± 0.0023(th), 
a = 0.5132 ± 0.0115(fit) ± 0.0381(th) , 

Compared to previous NLO results, the theoretical error on as{Mz) (which is dominated by the scale 
variation, improved at NNLO) is reduced by a factor of two, while the error on ao (which is dominated 
by the uncertainty on the Milan factor) remains unchanged. We observed that the sources of uncertainty 
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Figure 5: Perturbative and non-perturbative contributions to the moments of 1 — T at y/s — Mz as 
predicted by power corrections (left) and PYTHIA (right). 

on as(Mz) and ao largely decouple. An improvement on ao will only be achievable once the three- loop 
corrections to the Milan factor become available. 

It is noteworthy that application of the dispersive model to hadronization corrections results in a 
considerably lower value of a s (Mz) from event shapes than pervious studies based on Monte Carlo 
hadronization models [13], and in better agreement with measurements from other observables [59]. A 
direct comparison hints to an underestimation of hadronization effects in the Monte Carlo models. This 
feature has been observed previously also on the thrust distribution [47]. Revisiting the hadronization 
models in multi-purpose Monte Carlo programs appears to be mandatory for meaningful precision QCD 
studies at colliders. 
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A Tables of results 

In this appendix, we collect the extractions of a s (Mz) and ao at NLO and NNLO from individual 
moments of the six event shape variables: r, C, p, Y3, Bt, By/ . 



Pow. Corr.: ■ pert. non-pert 
Pythia: ■ pert. non-pert 
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NLO 


( T ) 


/-^2\ 
\ T / 


<- 3 ) 


<- 4 ) 




X Vdof 


1.0043 


1.0565 


0.8399 


0.6459 


0.4740 


otsiMz) 


0.1242 


0.1344 


0.1416 


0.1479 


0.1534 


Experimental Error 


0.0018 


0.0023 


0.0028 


0.0033 


0.0038 


variation: = 0.5 


-0.0054 


-0.0083 


-0.0101 


-0.0115 


-0.0129 


x„ = 2.0 


0.0066 


0.0102 


0.0123 


0.0143 


0.0162 


[ii variation: fij = 1.0 Gev 


0.0025 


0.0035 


0.0038 


0.0045 


0.0054 


= 3.0 Gev 


-0.0019 


-0.0025 


-0.0027 


-0.0031 


-0.0037 


M variation: M — 20 % 


0.0012 


0.0016 


0.0017 


0.0020 


0.0024 


.M + 20 % 


-0.0011 


-0.0014 


-0.0016 


-0.0018 


-0.0021 


Theoretical Error 


0.0072 


0.0109 


0.0130 


0.0151 


0.0172 




0.4782 


0.5147 


0.5359 


0.5521 


0.5744 


Experimental Error 


0.0151 


0.0152 


0.0189 


0.0222 


0.0243 


variation: x^ — 0.5 


0.0017 


0.0004 


-0.0038 


-0.0065 


-0.0081 


x„ = 2.0 


-0.0000 


-0.0001 


0.0030 


0.0051 


0.0064 


variation: M — 20 % 


0.0432 


0.0423 


0.0405 


0.0377 


0.0375 


A"! + 20 % 


-0.0306 


-0.0307 


-0.0298 


-0.0284 


-0.0290 


Theoretical Error 


0.0433 


0.0423 


0.0406 


0.0382 


0.0384 



Table 4: Results for a s (Q) and cto as obtained from fits to r moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NLO predictions. 



NNLO 


(r) 


(- 2 ) 


<- 3 > 


<- 4 ) 


(r 5 ) 


XVdof 


0.9889 


0.9411 


0.7284 


0.5526 


0.3997 


a s (M z ) 


0.1166 


0.1202 


0.1233 


0.1267 


0.1294 


Experimental Error 


0.0015 


0.0018 


0.0021 


0.0024 


0.0027 


x^ variation: x^ — 0.5 


-0.0020 


-0.0034 


-0.0042 


-0.0048 


-0.0054 


x^ = 2.0 


0.0025 


0.0042 


0.0051 


0.0058 


0.0065 


fix variation: fix = 1.0 Gev 


0.0017 


0.0017 


0.0017 


0.0019 


0.0022 


/Uj = 3.0 Gev 


-0.0011 


-0.0011 


-0.0011 


-0.0013 


-0.0014 


M variation: M - 20 % 


0.0009 


0.0010 


0.0009 


0.0011 


0.0012 


M + 20 % 


-0.0009 


-0.0009 


-0.0009 


-0.0010 


-0.0011 


Theoretical Error 


0.0032 


0.0046 


0.0054 


0.0062 


0.0070 


a 


0.5165 


0.5408 


0.5452 


0.5512 


0.5641 


Experimental Error 


0.0135 


0.0152 


0.0194 


0.0223 


0.0246 


x^ variation: x^ — 0.5 


0.0140 


0.0075 


0.0016 


-0.0008 


-0.0023 


a:,. = 2.0 


-0.0078 


-0.0045 


-0.0001 


0.0019 


0.0031 


M variation: M — 20 % 


0.0415 


0.0430 


0.0396 


0.0357 


0.0347 


M + 20 % 


-0.0298 


-0.0308 


-0.0286 


-0.0264 


-0.0261 


Theoretical Error 


0.0438 


0.0436 


0.0397 


0.0358 


0.0348 



Table 5: Results for a s (Q) and ao as obtained from fits to r moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NNLO predictions. 
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NLO 


(C) 


(c 2 ) 


(c 3 ) 


(° 4 ) 




X Vdof 


1.1849 


1.5245 


1.5651 


1.5446 


1.4094 


otsiMz) 


0.1230 


0.1308 


0.1347 


0.1374 


0.1407 


Experimental Error 


0.0013 


0.0016 


0.0020 


0.0023 


0.0026 


variation: = 0.5 


-0.0052 


-0.0079 


-0.0091 


-0.0100 


-0.0108 


x„ = 2.0 


0.0063 


0.0096 


0.0111 


0.0122 


0.0134 


[ii variation: fij = 1.0 Gev 


0.0029 


0.0045 


0.0051 


0.0057 


0.0064 


= 3.0 Gev 


-0.0022 


-0.0031 


-0.0034 


-0.0038 


-0.0041 


M variation: M — 20 % 


0.0013 


0.0020 


0.0022 


0.0025 


0.0028 


M + 20 % 


-0.0012 


-0.0018 


-0.0019 


-0.0022 


-0.0024 


Theoretical Error 


0.0071 


0.0107 


0.0124 


0.0137 


0.0151 




0.4267 


0.4632 


0.4789 


0.4839 


0.4857 


Experimental Error 


0.0082 


0.0064 


0.0067 


0.0069 


0.0070 


x^ variation: x^ — 0.5 


0.0052 


0.0027 


-0.0010 


-0.0035 


-0.0054 


x„ = 2.0 


-0.0029 


-0.0021 


0.0007 


0.0027 


0.0042 


variation: M — 20 % 


0.0324 


0.0359 


0.0377 


0.0376 


0.0366 


A"! + 20 % 


-0.0236 


-0.0266 


-0.0268 


-0.0283 


-0.0278 


Theoretical Error 


0.0328 


0.0360 


0.0377 


0.0377 


0.0370 



Table 6: Results for a s {Q) and «o as obtained from fits to C moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NLO predictions. 



NNLO 


(C) 


(c 2 ) 








XVdof 


1.1574 


1.2418 


1.2353 


1.1735 


1.0216 


a s (M z ) 


0.1161 


0.1180 


0.1193 


0.1202 


0.1216 


Experimental Error 


0.0011 


0.0013 


0.0016 


0.0017 


0.0019 


x^ variation: x^ — 0.5 


-0.0020 


-0.0033 


-0.0039 


-0.0043 


-0.0046 


x„ = 2.0 


0.0025 


0.0040 


0.0047 


0.0051 


0.0056 


fix variation: fix = 1.0 Gev 


0.0019 


0.0022 


0.0023 


0.0025 


0.0028 


/Uj = 3.0 Gev 


-0.0013 


-0.0014 


-0.0015 


-0.0016 


-0.0017 


M variation: M - 20 % 


0.0011 


0.0012 


0.0013 


0.0014 


0.0015 


M +20% 


-0.0010 


-0.0011 


-0.0012 


-0.0013 


-0.0014 


Theoretical Error 


0.0033 


0.0047 


0.0054 


0.0059 


0.0064 


a 


0.4689 


0.4897 


0.4919 


0.4877 


0.4828 


Experimental Error 


0.0071 


0.0063 


0.0067 


0.0069 


0.0070 


x^ variation: x^ — 0.5 


0.0166 


0.0095 


0.0053 


0.0027 


0.0010 


a:,. = 2.0 


-0.0105 


-0.0066 


-0.0033 


-0.0013 


0.0001 


M variation: M — 20 % 


0.0316 


0.0360 


0.0359 


0.0346 


0.0326 


M + 20 % 


-0.0234 


-0.0264 


-0.0265 


-0.0258 


-0.0246 


Theoretical Error 


0.0357 


0.0372 


0.0363 


0.0347 


0.0326 



Table 7: Results for a s (Q) and ao as obtained from fits to C moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NNLO predictions. 
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NLO 


(P) (p 2 ) 


x7dof 


0.6587 


0.7547 


a s (M z ) 


0.1164 


0.1152 


Experimental Error 


0.0023 


0.0033 


x u variation: x u — 0.5 


-0.0028 


-0.0038 


x a = 2.0 


0.0039 


0.0049 


pi variation: fij = 1.0 Gev 


0.0014 


0.0014 


Pi = 3.0 Gev 


-0.0011 


-0.0011 


M variation: M - 20 % 


0.0006 


0.0006 


M + 20 % 


-0.0006 


-0.0006 


Theoretical Error 


0.0042 


0.0051 


a 


0.5914 


0.5657 


Experimental Error 


0.0268 


0.0361 


a;^ variation: x M = 0.5 


0.0115 


0.0092 


x M - 2.0 


-0.0042 


-0.0047 


M variation: M — 20 % 


0.0795 


0.0748 


M + 20% 


-0.0539 


-0.0508 


Theoretical Error 


0.0803 


0.0753 



Table 8: Results for a s (Q) and «o as obtained from fits to p moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NLO predictions. 



NNLO 


(P) 


(P 2 ) 


xVdof 


0.6750 


0.7607 


a s (M z ) 


0.1142 


0.1113 


Experimental Error 


0.0021 


0.0030 


x^ variation: x M = 0.5 


-0.0009 


-0.0012 


Xfl = 2.0 


0.0013 


0.0017 


\ii variation: pj = 1.0 Gev 


0.0012 


0.0010 


pi = 3.0 Gev 


-0.0008 


-0.0007 


M variation: M - 20 % 


0.0007 


0.0006 


M + 20 % 


-0.0006 


-0.0006 


Theoretical Error 


0.0018 


0.0020 




0.6565 


0.6208 


Experimental Error 


0.0224 


0.0316 


x^ variation: x^ — 0.5 


0.0312 


0.0233 


x„ = 2.0 


-0.0184 


-0.0143 


M variation: M — 20 % 


0.0799 


0.0759 


M + 20 % 


-0.0547 


-0.0518 


Theoretical Error 


0.0858 


0.0794 



Table 9: Results for a s (Q) and ao as obtained from fits to p moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NNLO predictions. 
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NLO 


W (Yi) (Yi) (Yi) (Y*) 


rVdof 


0.8616 0.6386 0.7771 0.8691 0.9499 


a s {M z ) 


0.1183 0.1172 0.1165 0.1149 0.1124 


Experimental Error 


0.0011 0.0016 0.0020 0.0026 0.0033 


variation: x^ = 0.5 
x„ = 2.0 


-0.0040 -0.0042 -0.0040 -0.0038 -0.0035 
0.0053 0.0054 0.0052 0.0049 0.0045 


Theoretical Error 


0.0047 0.0048 0.0046 0.0043 0.0040 



Table 10: Results for a s (Q) and a as obtained from fits to Y 3 moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NLO predictions. 



NNLO 


(Y 3 ) 


(yi) 








x7dof 


0.8577 


0.6581 


0.7948 


0.8781 


0.9557 


a s {M z ) 


0.1156 


0.1136 


0.1136 


0.1129 


0.1106 


Experimental Error 


0.0010 


0.0015 


0.0019 


0.0025 


0.0032 


variation: x M = 0.5 


-0.0005 


-0.0008 


-0.0006 


-0.0002 


-0.0002 


x„ = 2.0 


0.0015 


0.0017 


0.0015 


0.0013 


0.0012 


Theoretical Error 


0.0010 


0.0013 


0.0011 


0.0008 


0.0007 



Table 11: Results for a s (Q) and ao as obtained from fits to I3 moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NNLO predictions. 



NLO 


(Bt) 




(53) 






X 2 /dof 


1.5775 


1.6741 


1.5926 


1.4005 


1.1996 


a s {M z ) 


0.1199 


0.1276 


0.1308 


0.1327 


0.1347 


Experimental Error 


0.0012 


0.0018 


0.0023 


0.0027 


0.0031 


x^ variation: x^ — 0.5 


-0.0037 


-0.0078 


-0.0093 


-0.0101 


-0.0108 


x„ = 2.0 


0.0049 


0.0094 


0.0112 


0.0123 


0.0133 


Hi variation: /ij = 1.0 Gev 


0.0021 


0.0032 


0.0036 


0.0038 


0.0041 


/Uj = 3.0 Gev 


-0.0016 


-0.0024 


-0.0026 


-0.0028 


-0.0030 


M variation: M - 20 % 


0.0010 


0.0015 


0.0016 


0.0017 


0.0018 


M + 20 % 


-0.0009 


-0.0014 


-0.0015 


-0.0016 


-0.0017 


Theoretical Error 


0.0054 


0.0101 


0.0119 


0.0130 


0.0140 


ao 


0.4252 


0.4897 


0.5180 


0.5193 


0.5088 


Experimental Error 


0.0130 


0.0105 


0.0112 


0.0129 


0.0146 


x^ variation: x^ — 0.5 


0.0154 


0.0083 


0.0005 


-0.0050 


-0.0093 


x» = 2.0 


-0.0106 


-0.0074 


-0.0015 


0.0031 


0.0070 


M variation: M — 20 % 


0.0333 


0.0443 


0.0500 


0.0493 


0.0452 


M+2G% 


-0.0238 


-0.0318 


-0.0358 


-0.0354 


-0.0329 


Theoretical Error 


0.0367 


0.0451 


0.0501 


0.0496 


0.0462 



Table 12: Results for a s {Q) and ao as obtained from fits to Bt moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NLO predictions. 



17 



NNLO 


(B T ) 


( B t) 


(B 3 T ) 


(B T ) 


( B t) 


X Vdof 


1.6191 


1.4765 


1.3723 


1.2059 


1.0363 


otsiMz) 


0.1164 


0.1158 


0.1158 


0.1158 


0.1162 


Experimental Error 


0.0011 


0.0014 


0.0017 


0.0019 


0.0022 


variation: = 0.5 


-0.0012 


-0.0027 


-0.0033 


-0.0035 


-0.0037 


x„ = 2.0 


0.0016 


0.0035 


0.0041 


0.0044 


0.0047 


[ii variation: fij = 1.0 Gev 


0.0017 


0.0017 


0.0018 


0.0018 


0.0019 


= 3.0 Gev 


-0.0012 


-0.0011 


-0.0012 


-0.0012 


-0.0012 


M variation: M — 20 % 


0.0010 


0.0010 


0.0010 


0.0010 


0.0011 


M + 20 % 


-0.0009 


-0.0009 


-0.0009 


-0.0009 


-0.0010 


Theoretical Error 


0.0025 


0.0040 


0.0046 


0.0049 


0.0052 




0.4844 


0.5053 


0.5059 


0.4938 


0.4772 


Experimental Error 


0.0104 


0.0094 


0.0098 


0.0108 


0.0117 


variation: x^ — 0.5 


0.0491 


0.0272 


0.0190 


0.0142 


0.0109 


x„ = 2.0 


-0.0295 


-0.0186 


-0.0129 


-0.0093 


-0.0066 


variation: M — 20 % 


0.0325 


0.0419 


0.0436 


0.0415 


0.0374 


A"! + 20 % 


-0.0240 


-0.0300 


-0.0312 


-0.0297 


-0.0270 


Theoretical Error 


0.0589 


0.0500 


0.0476 


0.0438 


0.0390 



Table 13: Results for a s (Q) and ao as obtained from fits to Bt moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NNLO predictions. 



NLO 


(B w ) 


(B 2 W ) 


(B 3 W ) 




(B 5 W ) 


XVdof 


1.5082 


1.2870 


1.1182 


0.8965 


0.6999 


a s (M z ) 


0.1128 


0.1077 


0.1049 


0.1023 


0.1010 


Experimental 


Error 


0.0015 


0.0020 


0.0027 


0.0033 


0.0039 


x^ variation: 


X[i — 


0.5 


0.0007 


-0.0028 


-0.0026 


-0.0022 


-0.0019 




Xfi = 


2.0 


0.0006 


0.0036 


0.0035 


0.0030 


0.0027 


fix variation: 


M/ = 


1.0 Gev 


0.0018 


0.0015 


0.0016 


0.0016 


0.0016 




A*/ = 


3.0 Gev 


-0.0014 


-0.0012 


-0.0012 


-0.0012 


-0.0013 


M. variation: 


M - 


20% 


0.0008 


0.0007 


0.0007 


0.0007 


0.0008 




M + 


20% 


-0.0008 


-0.0007 


-0.0007 


-0.0007 


-0.0007 


Theoretical Error 


0.0021 


0.0040 


0.0039 


0.0035 


0.0032 




0.3960 


0.3552 


0.3090 


0.2550 


0.2025 


Experimental 


Error 


0.0106 


0.0132 


0.0154 


0.0180 


0.0203 


x^ variation: 


Xfj, 


0.5 


0.0870 


0.0256 


0.0137 


0.0107 


0.0089 




Xfi = 


2.0 


-0.0401 


-0.0166 


-0.0097 


-0.0074 


-0.0059 


M variation: 


M - 


20% 


0.0357 


0.0308 


0.0222 


0.0112 


-0.0005 




M + 


20% 


-0.0250 


-0.0214 


-0.0157 


-0.0084 


-0.0006 


Theoretical Error 


0.0941 


0.0400 


0.0261 


0.0155 


0.0089 



Table 14: Results for a s (Q) and ao as obtained from fits to Bw moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NLO predictions. 
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NNLO 


(%} 










X 2 /dof 


1.5645 


1.2884 


1.1158 


0.8943 


0.7145 


a s (M z ) 


0.1117 


0.1058 


0.1032 


0.1014 


0.1037 


Experimental Error 


0.0014 


0.0018 


0.0025 


0.0031 


0.0038 


Xfj, variation: = 0.5 


-0.0010 


-0.0007 


-0.0006 


-0.0003 


0.0005 


x„ = 2.0 


0.0007 


0.0011 


0.0010 


0.0007 


0.0001 


fix variation: fix = 1.0 Gev 


0.0015 


0.0012 


0.0012 


0.0012 


0.0016 


/Uj = 3.0 Gev 


-0.0010 


-0.0008 


-0.0008 


-0.0008 


-0.0011 


Al variation: M - 20 % 


0.0009 


0.0007 


0.0007 


0.0007 


0.0009 


.M + 20 % 


-0.0008 


-0.0006 


-0.0007 


-0.0007 


-0.0009 


Theoretical Error 


0.0020 


0.0018 


0.0017 


0.0016 


0.0020 




0.4632 


0.4029 


0.3519 


0.2992 


0.2744 


Experimental Error 


0.0083 


0.0110 


0.0127 


0.0143 


0.0137 


rr^ variation: = 0.5 


0.0393 


0.0313 


0.0256 


0.0221 


0.0358 


^ = 2.0 


-0.0384 


-0.0232 


-0.0181 


-0.0155 


-0.0220 


M variation: M - 20 % 


0.0378 


0.0314 


0.0227 


0.0122 


0.0045 


X + 20 % 


-0.0270 


-0.0222 


-0.0163 


-0.0093 


-0.0045 


Theoretical Error 


0.0545 


0.0443 


0.0342 


0.0252 


0.0360 



Table 15: Results for a s (Q) and ao as obtained from fits to Bw moments measured by JADE and OPAL 
for centre-of-mass energies between 14.0 and 206.6 GeV using theoretical NNLO predictions. 
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